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Abstract. Matter-wave interferometry has been largely studied in the last few years. 
Usually, the main problem in the analysis of the diffraction experiments is to establish 
the causes for the loss of coherence observed in the interference pattern. In this 
work, we use different type of environmental couplings to model a two slit diffraction 
experiment with massive particles. For each model, we study the effects of decoherence 
on the interference pattern and define a visibility function that measures the loss of 
contrast of the interference fringes on a distant screen. Finally, we apply our results 
to the experimental reported data on massive particles C 70 . 
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The wave-particle duality of material objects is a hallmark of quantum mechanics. 
Up to now, the wave nature of particles has been demonstrated for electrons, neutrons, 
^ ■ atoms, and coherent atomic ensembles. Yet more important, many theoretical studies 
have been done around the mesoscopic systems [TJ [2]. Mesoscopic objects are neither 
microscopic nor macroscopic. They are generally systems that can be described by a 
^ ■ wavefunction, yet they are made up of a significant number of elementary constituents, 
such as atoms. Well-known examples these days are fullerene molecules Cqq and 
C70, which are expected to behave like classical particles. Nonetheless, the quantum 
interference of these molecules has been observed [3]. In particular, the advances with 
large molecules have stimulated the question what determines the limits to observe 
quantum delocalization with massive objects. Thus, there is a need to theoretically 
quantify the effect of decoherence (or dephasing) on the observed interference pattern 
in a double-slit experiment. It is quite intuitive that the resulting pattern shall be an 
interplay between the strength of the coupling to the environment, the slit separation 
and the distance the particle travels from the slit to the screen. 

We shall study a two slit diffraction experiment with particles of mass M diffracted 
by a grating (in the x direction) and then detected on a distant screen (located a 
distance L in the y direction). Note that coherence in the x direction is required in 
order to observe an interference pattern on the screen, whereas the dynamics in the y 
direction can be that of a free non- interacting particle since it just serves to transport 
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the particle from the grating to the screen. Initially, we may reasonably assume that the 
action of the grating is to prepare a superposition of two Gaussian wavepackets (which 
best describe a massive particle), centered at each location of the respective slits and 
factorized as [HE] 0) = (</> x (x, 0) + (/) 2 (x, 0))x(y, 0) <g> ((X, 0), where |0i| 2 and |0 2 | 2 
correspond to the probability amplitudes for the particle to pass through slit 1 and slit 
2 (in the x-axis), respectively, while xiv^t) represents the Gaussian wave function in 
the y direction (where no superposition is needed) and ((X,t) describes the state of 
the environment coupled to the subsystem. Note that we are assuming translational 
invariance in the z-axis [5]. 

We shall consider the environment as a set of non-interacting harmonic oscillators 
and the dynamics of the particles modeled by a quantum brownian motion (QBM) [6j. 
This behavior can be imputed to the passing of the particles through the slits producing 
vibrations or any other kind of interactions with the walls of the grating able to corrupt 
the interference pattern. Moreover, also the finite size of the grating and the differences 
in the slit apertures can attenuate the visibility of the interference fringes, especially in 
the case of complex molecules such us C 60 and C 70 [7\ . In order to study the interference 
pattern registered on the screen at a later time ti-, we need to obtain the evolution in 
time of the reduced density matrix p r (x,x',t), which is given by the following master 
equation 

dp r (x,x',t) ih fd 2 p r d 2 p r \ T>(t) , 2 

dt = 2M\dx^ ~ dx 1 *) ~^f^ x ~ x) Pr 

where 7(i) is the dissipative coefficient (proportional to the square of the coupling 
constant to the environment), T>(t) the diffusive coefficient and f(t) the coefficient 
responsible for the anomalous diffusion. Eq.([T]) has been obtained by assuming the 
environment to be in equilibrium, at a temperature T. In the case that the system is 
coupled to an ohmic environment in the high temperature limit (fc^T >> hu), these 
coefficients are constant j(t) = 70, T>(t) = 2M r y kBT and f(t) ~ 1/ksT in units of 
h = 1 [BJ. It is important to stress that the equation of movement for the generally used 
scattering model i^f- = [H, p r ] — iA[x, [x, p r \], where the effect of the environment is 
included in the collision term A, can be obtained from Eq.([T]) in the high temperature 
limit (neglecting dissipation) for the markovian case if written in the Lindblad form. 
However, this master equation refers to a more general movement that can be used for all 
temperatures, even to study the dynamics of the test particle at zero temperature. Yet 
more interesting, this formulation includes the phenomenological scattering model and 
verifies the fluctuation-dissipation theorem for a general system in thermal equilibrium. 

The interference pattern corresponds to the probability distribution of the time 
evolved wave function: 

P(x,t) = \^(x,t)\ 2 = L 1 (x,t)*(i) 1 {x,t)+<i) 2 {x,ty<f) 2 {x,t) 

+ <f) 2 (x,t)*(f> 1 (x,t) + (t> 1 (x,t)*4> 2 (x,t)) ®x(y,t)*x(y',t), 
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Figure 1. (a) Evolution in time for the interference pattern registered on the screen 
when the system is coupled to an environment. Parameters used: Lq = 2, <j x q = 0.5, 
M = 1, 70 = 0.001 and fc^T = 300. (b) For the same values the wigner function at 
t = 2 s shows that the interferences have disappeared since the wigner is positive. 

which is proportional to the diagonal part of the density matrix defined as p(x, x', t) = 
\^f(x,i))(^f(x',t)\. In the case we are studying herein, i.e an open quantum system, the 
interference pattern at a given time t on the screen is given by: 

P(x , t) = t)\ 2 + |0 2 (x, t)\ 2 + 2T(t)Re( ( f)* 1 (x, t)(j) 2 (x, t))) \ X (y, t)\ 2 , 

where the terms in brackets correspont to p r (x,x,t) and T(t) encodes the information 
about the statistical nature of noise since it is obtained after tracing out the degrees of 
freedom of the environment. T(t) is an exponential decaying term which suppresses 
the interference terms in a decoherence time scale to- In the case of the QBM, 
T(t) = e~ Ax fo v ( s ) ds (with Ax 2 = (x — x') 2 ) represents the noise induced environmental 
effect on the system. In this framework, we numerically solved the master equation 
Eq.flU) for the reduced density matrix t) of two well localized Gaussian wave 

packets initially given by 

0) = iv(exp( (3 7f o)2 ) + exp( (x + f o)2 )) ® expi--^ - tk y y), 

where 2Lo is the initial separation of the center of the wave packets, a 2 and a 2 are the 
initial width of the packet in the x and y-axis, respectively, and k y the initial moment of 
the particle in the y direction. It is important to note that Lq, a x o, cr y o and k y are all free 
parameters that have to be tuned with the experimental data. In addition, we assume 
that Ap y << p y , so the moment component is sharply defined and the wave packet has 
a characteristic wavelength associated X^b ~ ft/Py « Ay [Hl[9]. Once the reduced 
density matrix is known for all time t, we know the dynamics of these two packets and 
can study the effects of decoherence on the interference pattern registered on the distant 
screen. In FigJT](a) we present the time evolution for the interference pattern registered 
on the screen. We can there observe the two wave packets initially separated a distance 
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Figure 2. (a) Time evolution for the visibility fringe v{t) for ksT — 300, 70 = 
0.001, Lq = 2, M = 1 and a x o = 0.5. The estimation of the decoherence time 
to ~ l/(2MjokBTLQ) — 0.41 s coincides with the timescale at which the visibility 
starts to decrease towards a null value, (b) Time evolution for the visibility function 
vcif) for neutrons (Cn =0.1) and fullerenes (Cf = 1 and Cf = 2 ) in the presence of 
an external classical time dependent electromagnetic field (incoherence effects). The 
curves are for the first minimum and maximum of the interference pattern and reach 
a different asymptotic value than in the case of v(t). 



2Lq which start to spread as time increases. In a short timescale, interferences start to 
develop, however we can clearly observe that the minimum fringes of the interference 
pattern are not exactly zero. This means that there is a loss of contrast between theis 
case (an open system) with respect to the interference pattern of an unitary evolution. 
This means that decoherence is already playing a crucial role in determining the wave 
behavior of the massive particle. In Figfjjb) we present the wigner function for the 
corresponding reduced density matrix at a fixed time t = 2 s. We can see that it is a 
completely positive function, which assures us that decoherence has been effective in this 
short timescale. We shall estimate the decoherence time tp, i.e. the timescale for which 
the interferences are mostly destroyed (up to a 70 %) as Tvitr,) = exp(—VAx 2 t) ~ 1/e. 
It is easily deduced that tp w 1/(2) Ax 2 ). A quantity of particular importance in matter 
wave interferometry is the fringe visibility v{t). In order to study the role of decoherence 
on the fringe visibility, we define it as v(t) ~ — \p in t(x,x,t)\ — _ — — r(t] — — where 

J ' V ' Pll(X,X,t)+p22(X,X,t) pii(x,X,t)+p22(X,X,t) 

Pa = \4>i{%, £)| 2 j with i = 1,2 and p irA the interference terms. 

Clearly, the visibility fringe goes down as t^, i.e. the observation time, is larger than 
the decoherence time t^,. However, if we succeed in performing our two slit experiment 
in a time ti < to at a fixed room temperature ksT, we can see that the visibility 
fringe decreases as the coupling to the environment (70) increases. This is so, because 
the decoherence time depends inversely on the coupling constant: the bigger 70, the 
shorter the decoherence time. Not only can we check the dependence upon the coupling 
constant but also its time evolution. This behavior is shown in Figj2]^a). For short 
times, the visibility increases from zero to a maximum value because the interferences 
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Figure 3. The interference pattern (y ~ 0.62) registered on the screen considering 
the incoherence (or dephasing) effects to model a two slit diffraction experiment with 
massive particles Cro- The curves are for the isolated case, the theoretical prediction 
using the experimental data reported in [3] and the incoherence model using Tc with 
(7 = 1. 

start to develop at that short timescale but are not present at t — (since the wave 
packets are initially separated and have to spread so as to generate the interferences 
see FigJT](a)). This maximum value coincides with the estimated decoherence time tp. 
Then, the visibility starts to decrease, since the destruction of the interferences is taking 
place. Note that the visibility is a quantity that measures the loss of contrast of the 
interference fringes. Then, it is expected that those with the bigger contrast suffer from 
this attenuation the more, as seen in FigfS^a). Clearly, the observation time t^ must be 
shorter than the decoherence time in order to observe the interference pattern. 

It shall also be interesting to study the visibility function for a model of incoherence 
such as the one presented by us in [10]. In particular, T c = Jq(\C\) is constant in time for 
the experimental data of neutrons and fullerenes pU] ■ Therein, we previously estimated 
the quantity C for these massive particles and observed that, contrary to might be 
naively expected, in thought and real experiments, Cf ~ 0(1). However, for neutral 
particles with permanent dipole moment this value is much lower Cn ~ (9(0.01)— 0(0.1). 
Therefore, we define the visibility function vn(t) as vAt) = — -, — J °(\ c ^ — - and show 
its time evolution for neutrons and fullerenes in Figj2](b). We recall that this interaction 
is always present in charged and neutral particles with permanent dipole moment and 
can never be turned off. However, in the case of neutrons it has been shown in (TO] that 
it can be neglected, whereas it is unexpectedly important for massive particles such as 
fullerenes. 

Finally, we shall apply our models of decoherence and incoherence to the 
experimental data reported in [3] to reproduce the observed pattern for fullerenes CV . 
In FigBJ the interference pattern is shown. Therein, we have considered the unitary and 
nonunitary evolution of the particles. For these massive particles, we can see that the 
interference pattern is attenuated when the system is open. What is more important, is 
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that the modeling of the incoherence effects through the overlap factor Tc is enough to 
reproduce the effects of the environment on the interference pattern of a real diffraction 
experiment with massive particles. 

All in all, we have presented a fully quantum mechanical treatment using 
a microscopic model of environment and also a concrete example to include the 
incoherence effects. Therefore, we have studied the effects of decoherence on the 
interference pattern of thought diffraction experiments and presented an analysis of 
matter wave interferometry in the presence of a dynamical quantum environment 
such as the quantum brownian motion model. We have shown the interference 
patterns and visibility function v(t) for thought diffracted free particles in the 
high temperature limit (assumption valid for massive particles diffracted at room 
temperature). Yet more important, we have defined the visibility function vc{t) for 
a model environment previously developed which describes the incoherence effects 
originated in the experimental difficulty of producing the same initial/final state for all 
particles (i.e the existence of a random variable such as the particle's emission time). We 
showed that it is qualitatively different than the one commonly found in the literature 
and very important in the case of diffraction experiments with massive particles such 
as fullerenes. In this case, the incoherence effects are enough to model the attenuation 
of the interference pattern observed in the real experiment, whereas in the case of cold 
neutrons the incoherence effects are not of such importance [TT]. Therefore, in the latter 
case we must consider the decoherence effects by using the corresponding formulation 
(such us QBM). This result might have been expected since the interaction of more 
massive particles with the external classical field is more important than for those with 
a smaller mass where other kind of interactions (such as with the walls of the grating 
or the air molecules) seem to prevail. 
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